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Abstract 

Recent experiments on conduction between a semiconductor and a super- 
conductor have revealed a variety of new mesoscopic phenomena. Here is 
a review of the present status of this rapidly developing field. A scattering 
theory is described which leads to a conductance formula analogous to Lan- 
dauer's formula in normal-state conduction. The theory is used to identify 
features in the conductance which can serve as "signatures" of phase-coherent 
Andreev reflection, i.e. for which the phase coherence of the electrons and the 
Andreev-reflected holes is essential. The applications of the theory include 
a quantum point contact (conductance quantization in multiples of 4e 2 //i), 
a quantum dot (non-Lorentzian conductance resonance), and quantum in- 
terference effects in a disordered normal-superconductor junction (enhanced 
weak-localization and reflectionless tunneling through a potential barrier). 
The final two sections deal with the effects of Andreev reflection on universal 
conductance fluctuations and on the shot noise. 
Lectures at the Les Houches summer school, Session LXI, 
1994, to be published in: Mesoscopic Quantum Physics, E. Akkermans, 
G. Montambaux, and J.-L. Pichard, eds. (North-Holland, Amsterdam). 
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FIG. 1. Normal reflection by an insulator (I) versus Andreev reflection by a superconductor (S) 
of an electron excitation in a normal metal (N) near the Fermi level. Normal reflection (left) con- 
serves charge but does not conserve momentum. Andreev reflection (right) conserves momentum 
but does not conserve charge: The electron (e) is reflected as a hole (h) with the same momentum 
and opposite velocity. The missing charge of 2e is absorbed as a Cooper pair by the superconducting 
condensate. 



At the interface between a normal metal and a superconductor, dissipative electrical 
current is converted into dissipationless supercurrent. The mechanism for this conversion 
was discovered thirty years ago by A. F. Andreev [1]: An electron excitation slightly above 
the Fermi level in the normal metal is reflected at the interface as a hole excitation slightly 
below the Fermi level (see fig. 1). The missing charge of 2e is removed as a supercurrent. 
The reflected hole has (approximately) the same momentum as the incident electron. (The 
two momenta are precisely equal at the Fermi level.) The velocity of the hole is minus the 
velocity of the electron (cf. the notion of a hole as a "time-reversed" electron). This curious 
scattering process is known as retro-reflection or Andreev reflection. 

The early theoretical work on the conductance of a normal-metal - superconductor (NS) 
junction treats the dynamics of the quasiparticle excitations semiclassically, as is appropri- 
ate for macroscopic junctions. Phase coherence of the electrons and the Andreev-reflected 
holes is ignored. Interest in "mesoscopic" NS junctions, where phase coherence plays an 
important role, is a recent development. Significant advances have been made during the 
last few years in our understanding of quantum interference effects due to phase-coherent 
Andreev reflection. Much of the motivation has come from the technological advances in 
the fabrication of a highly transparent contact between a superconducting film and the two- 
dimensional electron gas in a semiconductor heterostructure. These systems are ideal for 
the study of the interplay of Andreev reflection and the mesoscopic effects known to occur 
in semiconductor nanostructures [2] , because of the large Fermi wavelength, large mean free 
path, and because of the possibility to confine the carriers electrostatically by means of gate 
electrodes. In this series of lectures we review the present status of this rapidly developing 
field of research. 

To appreciate the importance of phase coherence in NS junctions, consider the resistance 
of a normal-metal wire (length L, mean free path /). This resistance increases monotonically 
with L. Now attach the wire to a superconductor via a tunnel barrier (transmission probabil- 
ity T). Then the resistance has a minimum when L ~ l/Y. The minimum disappears if the 
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phase coherence between the electrons and holes is destroyed, by increasing the voltage or 
by applying a magnetic field. The resistance minimum is associated with the crossover from 
a r _1 to a T~ 2 dependence on the barrier transparency. The T~ 2 dependence is as expected 
for tunneling into a superconductor, being a two-particle process. The T _1 dependence is 
surprising. It is as if the Andreev-reflected hole can tunnel through the barrier without 
reflections. This socalled "reflectionless tunneling" requires relatively transparent NS inter- 
faces, with T > l/L. Semiconductor — superconductor junctions are convenient, since the 
Schottky barrier at the interface is much more transparent than a typical dielectric tunnel 
barrier. The technological effort is directed towards making the interface as transparent as 
possible. A nearly ideal NS interface (T ~ 1) is required if one wishes to study how Andreev 
reflection modifies the quantum interference effects in the normal state. (For T 1 these 
are obscured by the much larger reflectionless-tunneling effect.) The modifications can be 
quite remarkable. We discuss two examples. 

The first is weak localization. In the normal state, weak localization can not be detected 
in the current-voltage (I-V) characteristic, but requires application of a magnetic field. The 
reason is that application of a voltage (in contrast to a magnetic field) does not break time- 
reversal symmetry. In an NS junction, however, weak localization can be detected in the 
I-V characteristic, because application of a voltage destroys the phase coherence between 
electrons and holes. The result is a small dip in dl/dV versus V around V — for F ~ 1. 
On reducing T, the dip crosses over to a peak due to reflectionless tunneling. The peak is 
much larger than the dip, but the widths are approximately the same. 

The second example is universal conductance fluctuations. In the normal state, the 
conductance fluctuates from sample to sample with a variance which is independent of 
sample size or degree of disorder. This is one aspect of the universality. The other aspect 
is that breaking of time-reversal symmetry (by a magnetic field) reduces the variance by 
precisely a factor of two. In an NS junction, the conductance fluctuations are also size and 
disorder independent. However, application of a time- reversal-symmetry breaking magnetic 
field has no effect on the magnitude. 

These three phenomena, weak localization, reflectionless tunneling, and universal con- 
ductance fluctuations, are discussed in sections 4, 5, and 6, respectively. Sections 2 and 3 
are devoted to a description of the theoretical method and to a few illustrative applications. 
The method is a scattering theory, which relates the conductance Gns of the NS junction to 
the N x N transmission matrix t in the normal state (N is the number of transverse modes 
at the Fermi level). In the limit of zero temperature, zero voltage, and zero magnetic field, 
the relationship is 



where the transmission eigenvalue T n is an eigenvalue of the matrix product ttf . The same 
numbers T n (n — 1, 2, . . . N) determine the conductance Gn in the normal state, according 
to the Landauer formula 




(1.1) 




2 N 



(1.2) 
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The fact that the same transmission eigenvalues determine both and Gns means that one 
can use the same (numerical and analytical) techniques developed for quantum transport in 
the normal state. This is a substantial technical and conceptual simplification. 

The scattering theory can also be used for other transport properties, other than the 
conductance, both in the normal and the superconducting state. An example, discussed in 
section 7, is the shot noise due to the discreteness of the carriers. A doubling of the ratio 
of shot-noise power to current can occur in an NS junction, consistent with the notion of 
Cooper pair transport in the superconductor. 

We conclude in section 8. 

We restrict ourselves in this review (with one exception) to two-terminal geometries, 
with a single NS interface. Equation (1.1), as well as the Landauer formula (1.2), only 
describe the two-terminal conductance. More complex multi-terminal geometries, involving 
several NS interfaces, have been studied theoretically by Lambert and coworkers [3,4], and 
experimentally by Petrashov et al. [5]. Since we focus on phase-coherent effects, most of 
our discussion concerns the linear-response regime of infinitesimal applied voltage. A recent 
review by Klapwijk contains a more extensive coverage of the non-linear response at higher 
voltages [6]. The scattering approach has also been applied to the Josephson effect in SNS 
junctions [7] , resulting in a formula for the super current-phase relationship in terms of the 
transmission eigenvalues T n in the normal state. We do not discuss the Josephson effect 
here, but refer to ref. [8] for a review of mesoscopic SNS junctions. Taken together, ref. [8] 
and the present work describe a unified approach to mesoscopic superconductivity. 

II. SCATTERING THEORY 

The model considered is illustrated in fig. 2. It consists of a disordered normal region 
(hatched) adjacent to a superconductor (S). The disordered region may also contain a ge- 
ometrical constriction or a tunnel barrier. To obtain a well-defined scattering problem we 
insert ideal (impurity-free) normal leads N x and N 2 to the left and right of the disordered 
region. The NS interface is located at x — 0. We assume that the only scattering in the 
superconductor consists of Andreev reflection at the NS interface, i.e. we consider the case 
that the disorder is contained entirely within the normal region. The spatial separation of 
Andreev and normal scattering is the key simplification which allows us to relate the con- 
ductance directly to the normal-state scattering matrix. The model is directly applicable 
to a superconductor in the clean limit (mean free path in S large compared to the super- 
conducting coherence length £), or to a point-contact junction (formed by a constriction 
which is narrow compared to £). In both cases the contribution of scattering within the 
superconductor to the junction resistance can be neglected [9]. 

The scattering states at energy e are eigenfunctions of the Bogoliubov-de Gennes (BdG) 
equation. This equation has the form of two Schrodinger equations for electron and hole 
wavefunctions u(f) and v(r), coupled by the pair potential A(r) [10]: 
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FIG. 2. Normal-metal-superconductor junction containing a disordered normal region 
(hatched). Scattering states in the two normal leads Ni and N2 are indicated schematically. 

Here H Q = (p + eA) 2 /2m + V - E F is the single-electron Hamiltonian, containing an electro- 
static potential V(r) and vector potential A(r). The excitation energy e is measured relative 
to the Fermi energy E F . To simplify construction of the scattering basis we assume that the 
magnetic field B (in the z-direction) vanishes outside the disordered region. One can then 
choose a gauge such that A = in lead N 2 and in S, while A X ,A Z = 0, A y — Ai = constant 
in lead 

The pair potential in the bulk of the superconductor (x ^> £) has amplitude A Q and 
phase 0. The spatial dependence of A(r) near the NS interface is determined by the self- 
consistency relation [10] 

m = l»(r1lE^(r)u(0[l-2/( e )], (2.2) 

£>0 

where the sum is over all states with positive eigenvalue, and f(e) = [1 + exp^/A^T)] -1 
is the Fermi function. The coefficient g is the interaction constant of the BCS theory of 
superconductivity. At an NS interface, g drops abruptly (over atomic distances) to zero, in 
the assumed absence of any pairing interaction in the normal region. Therefore, A(r) = 
for x < 0. At the superconducting side of the NS interface, A{r) recovers its bulk value 
Aoe 1<?i only at some distance from the interface. We will neglect the suppression of A(f) on 
approaching the NS interface, and use the step-function model 

A{f) = A o e i4, 0{x). (2.3) 

This model is also referred to in the literature as a "rigid boundary-condition". Likharev 
[11] discusses in detail the conditions for its validity: If the width W of the NS junction is 
small compared to £, the non-uniformities in A(f) extend only over a distance of order W 
from the junction (because of "geometrical dilution" of the influence of the narrow junction 
in the wide superconductor). Since non-uniformities on length scales <C £ do not affect the 
dynamics of the quasiparticles, these can be neglected and the step-function model holds. A 
point contact or microbridge belongs in general to this class of junctions. Alternatively, the 
step-function model holds also for a wide junction if the resistivity of the junction region is 
much bigger than the resistivity of the bulk superconductor. This condition is formulated 
more precisely in ref. [11]. A semiconductor — superconductor junction is typically in this 
second category. Note that both the two cases are consistent with our assumption that the 
disorder is contained entirely within the normal region. 

It is worth emphasizing that the absence of a pairing interaction in the normal region 
(g(r) = for x < 0) implies a vanishing pair potential A(r), according to eq. (2.2), but does 
not imply a vanishing order parameter \l/(r), which is given by 
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*(fO=Ev*(f)u(fO[l-2/( e )]. 

£>0 



(2.4) 



Phase coherence between the electron and hole wave functions u and v leads to ty(r) 7^ 
for x < 0. The term "proximity effect" can therefore mean two different things: One is the 
suppression of the pair potential A at the superconducting side of the NS interface. This 
is a small effect which is neglected in the present work (and in most other papers in this 
field) . The other is the induction of a non-zero order parameter \l> at the normal side of the 
NS interface. This effect is fully included here, even though \1/ does not appear explicitly in 
the expressions which follow. The reason is that the order parameter quantifies the degree 
of phase coherence between electrons and holes, but does not itself affect the dynamics of 
the quasiparticles. (The BdG equation (2.1) contains A not ^ .) 

We now construct a basis for the scattering matrix (s-matrix). In the normal lead N 2 
the eigenfunctions of the BdG equation (2.1) can be written in the form 

< e (N 2 ) = ( 1 )(K)- 1 / 2 K(y,z)exp(±iklx), 

1^h(N 2 ) = (?)(^)- 1/2 ^(y^)exp(±i^), (2.5) 

where the wavenumbers k* and k\ are given by 

kt = {2m/Ti 2 f 2 {E F - E n + a c ' h ef/\ (2.6) 

and we have defined a e = 1, cr h = — 1. The labels e and h indicate the electron or hole 
character of the wavefunction. The index n labels the modes, <& n (y,z) is the transverse 
wavefunction of the n-th mode, and E n its threshold energy: 

[(p 2 y + p 2 z )/2m + V(y, z)]<S> n (y, z) = E n <5> n (y, z). (2.7) 

The eigenfunction $ n is normalized to unity, Jdy fdz\§ n \ 2 = 1. With this normaliza- 
tion each wavefunction in the basis (2.5) carries the same amount of quasiparticle current. 
The eigenfunctions in lead Ni are chosen similarly, but with an additional phase factor 
exp[— ia c ' h (eAi/h)y] from the vector potential. 

A wave incident on the disordered normal region is described in the basis (2.5) by a 
vector of coefficients 

4° = (c e + (N 1 ),c e -(N 2 ),^(N 1 ),c+(N 2 )). (2.8) 

(The mode-index n has been suppressed for simplicity of notation.) The reflected and 
transmitted wave has vector of coefficients 

cT = (c e -(N 1 ),c e + (N 2 ),c+(N 1 ),c 1 ;(N 2 )). (2.9) 

The s-matrix s N of the normal region relates these two vectors, 

c° N ut = %4 n . (2.10) 
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Because the normal region does not couple electrons and holes, this matrix has the block- 
diagonal form 



s N (e) = 



so(e) 
s ( 








) 




^21 r 22 



)■ 



(2.11) 



Here so is the unitary s-matrix associated with the single-electron Hamiltonian H,q. The 
reflection and transmission matrices r(e) and t(e) are N x N matrices, N(e) being the 
number of propagating modes at energy e. (We assume for simplicity that the number of 
modes in leads Ni and N 2 is the same.) The matrix s is unitary (sqS = 1) and satisfies 
the symmetry relation s (e,B)ij = sq(e, —B)ji. 

For energies < e < A Q there are no propagating modes in the superconductor. We can 
then define an s-matrix for Andreev reflection at the NS interface which relates the vector 
of coefficients (c~(N 2 ), cjJ~(N 2 )) to (c+(N 2 ), (N 2 )). The elements of this s-matrix can be 
obtained by matching the wavefunctions (2.5) at x = to the decaying solutions in S of the 
BdG equation. If terms of order A /E F are neglected (the socalled Andreev approximation 
[1]), the result is simply 



where a = exp[— iarccos(e/^ )]. Andreev reflection transforms an electron mode into a hole 
mode, without change of mode index. The transformation is accompanied by a phase shift, 
which consists of two parts: 

1. A phase shift — arccos(e/Z\ ) due to the penetration of the wavefunction into the 
superconductor. 

2. A phase shift equal to plus or minus the phase of the pair potential in the supercon- 
ductor (plus for reflection from hole to electron, minus for the reverse process). 

We can combine the 2N linear relations (2.12) with the 4N relations (2.10) to obtain a 
set of 2N linear relations between the incident wave in lead N x and the reflected wave in the 
same lead: 



The four N x N matrices s ee , Shh, Sch, and She form together the scattering matrix s of the 
whole system for energies < e < A . An electron incident in lead N\ is reflected either as 
an electron (with scattering amplitudes s ee ) or as a hole (with scattering amplitudes Sh e )- 
Similarly, the matrices Shh and s c h contain the scattering amplitudes for reflection of a hole 
as a hole or as an electron. After some algebra we find for these matrices the expressions 



c e -(N 2 )=ae\-(N 2 ), 
c h + (N 2 ) = ae-\ + (N 2 ), 



(2.12) 



c e -(N 1 ) = s ee c+(N 1 ) + s eh Ch(N 1 ), 
c+(N 1 ) = s hc c c + (N 1 ) + s hh Ch(N 1 ). 



(2.13) 



See(^) = + Oc\ 2 (£)r^ (-£■ )M c t 21 (£T ), 

Shh(^) = rn(-e) + a 2 tl 2 (-e)r 22 (e)M h t* 21 (-e), 
s ch (e) = ae i % 2 (e)M h t* 21 (-e), 
s hc (e) = ae-' l Hl 2 (-e)M c t 2l (e), 



(2.14) 
(2.15) 
(2.16) 
(2.17) 
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where we have defined the matrices 



M e = [1 - a 2 r 22 ( £ )r 2 * 2 (-£)]-\ 

M h = [1 - a 2 r* 22 {-s)r 22 {s)\-\ (2.18) 

One can verify that the s-matrix constructed from these four sub-matrices satisfies uni- 
tarity (s^s = 1) and the symmetry relation s(e,B,(f))ij = s(e, — B, — cf))^, as required by 
quasiparticle-current conservation and by time-reversal invariance, respectively. 

For the linear-response conductance Gns of the NS junction at zero temperature we only 
need the s-matrix at the Fermi level, i.e. at e = 0. We restrict ourselves to this case and 
omit the argument e in what follows. We apply the general formula [12-14] 

2e 2 , 4 e 2 

Cns = -^-Tr (1 - s ce 4 + s h Xe) = — Tr s he s{ e - ( 2 - 19 ) 

The second equality follows from unitarity of s, which implies 1 — s ee sl e = s eh s\ h = 
(-4e)~ ls L s hc s L> so that Tr (1 — s ee sl e ) = Tr s he -sL- We now substitute eq. (2.17) for e = 
(a = — i) into eq. (2.19), and obtain the expression 

4e 2 

G NS = — Tr4t 12 (l + r* 22 r 22 )-H* 21 tl(l + r 22 r 2 T 2 )-\ (2.20) 

where M T = (M*)t denotes the transpose of a matrix. The advantage of eq. (2.20) over 
eq. (2.19) is that the former can be evaluated by using standard techniques developed for 
quantum transport in the normal state, since the only input is the normal-state scattering 
matrix. The effects of multiple Andreev reflections are fully incorporated by the two matrix 
inversions in eq. (2.20). 

In the absence of a magnetic field the general formula (2.20) simplifies considerably. 
Since the s-matrix sq of the normal region is symmetric for B = 0, one has r 22 = r 22 and 
t l2 = t 21 . Equation (2.20) then takes the form 

4e 2 

Cns = -^-Trt f 12 t 12 (l + 4 2 r 22 )- 1 t t 12 t 12 (l + r\ 2 r 22 Y x 

= ^Tr (4t 12 (2 - A^y 1 ) 2 . (2.21) 

In the second equality we have used the unitarity relation r 22 r 22 + t\ 2 t 12 = 1. The trace 
(2.21) depends only on the eigenvalues of the Hermitian matrix t[ 2 t l2 . We denote these 
eigenvalues by T n (n = 1, 2, ... N). Since the matrices t\ 2 t 12 , t l2 t{ 2 , t 21 t 21 , and t 21 t 2l all have 
the same set of eigenvalues, we can omit the indices and write simply ttf . We obtain the 
following relation between the conductance and the transmission eigenvalues: 

4e 2 N T 2 

Gns = xS(23%f- (2 - 22) 

This is the central result of ref. [15]. 

Equation (2.22) holds for an arbitrary transmission matrix t, i.e. for arbitrary disorder 
potential. It is the multi-channel generalization of a formula first obtained by Blonder, 
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Tinkham, and Klapwijk [12] (and subsequently by Shelankov [16] and by Zaitsev [17]) for 
the single- channel case (appropriate for a geometry such as a planar tunnel barrier, where 
the different scattering channels are uncoupled). A formula of similar generality for the 
normal-metal conductance G^ is the multi-channel Landauer formula 

2p 2 2p 2 n 

Gn = —Trtt* = — ]TT n . (2.23) 

11 11 71=1 

In contrast to the Landauer formula, eq. (2.22) for the conductance of an NS junction 
is a non-linear function of the transmission eigenvalues T n . When dealing with a non- 
linear multi-channel formula as eq. (2.22), it is of importance to distinguish between the 
transmission eigenvalue T n and the modal transmission probability T n = J2m=i \tnm\ 2 - The 
former is an eigenvalue of the matrix tt\ the latter a diagonal element of that matrix. The 
Landauer formula (2.23) can be written equivalently as a sum over eigenvalues or as sum 
over modal transmission probabilities: 

h N N 

-G N = ET^E^- (2-24) 

Ze 71=1 71=1 

This equivalence is of importance for (numerical) evaluations of the Landauer formula, in 
which one calculates the probability that an electron injected in mode n is transmitted, and 
then obtains the conductance by summing over all modes. The non-linear scattering formula 
(2.22), in contrast, can not be written in terms of modal transmission probabilities alone: 
The off-diagonal elements of ttf contribute to Gns i n an essential way. Previous attempts 
to generalize the one-dimensional Blonder-Tinkham-Klapwijk formula to more dimensions 
by summing over modal transmission probabilities (or, equivalently, by angular averaging) 
were not successful precisely because only the diagonal elements of ttf were considered. 

III. THREE SIMPLE APPLICATIONS 

To illustrate the power and generality of the scattering formula (2.22), we discuss in this 
section three simple applications to the ballistic, resonant-tunneling, and diffusive transport 
regimes [15]. 

A. Quantum point contact 

Consider first the case that the normal metal consists of a ballistic constriction with 
a normal-state conductance quantized at G^ = 2N e 2 /h (a quantum point contact). The 
integer N Q is the number of occupied one- dimensional subbands (per spin direction) in the 
constriction, or alternatively the number of transverse modes at the Fermi level which can 
propagate through the constriction. Note that N <C N. An "ideal" quantum point contact 
is characterized by a special set of transmission eigenvalues, which are equal to either zero 
or one [2]: 
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FIG. 3. Solid curve: Conductance Gns versus Fermi energy of a quantum point contact between 
a normal and a superconducting reservoir (shown schematically in the inset). The dotted curve is 
twice the conductance Gn for the case of two normal reservoirs [20]. The constriction is defined by 
the 2D saddle-point potential V(x, y) =Vq — ^muj 2 x 2 + \rriLo 2 y 2 , with uj y /u) x = 3; Gns is calculated 
from eq. (2.22), with T n = [1 + exp(-27r£ n /?kj, E )]~ 1 , e n = Ep — Vq — (n — ^)hu y . (From ref. [15].) 

T.-Wl 1 *???* (3.1) 



if iVo < n < N, 



where the eigenvalues have been ordered from large to small. We emphasize that eq. (3.1) 
does not imply that the transport through the constriction is adiabatic. In the case of 
adiabatic transport, the transmission eigenvalue T n is equal to the modal transmission prob- 
ability T n . In the absence of adiabaticity there is no direct relation between T n and T n . 
Substitution of eq. (3.1) into eq. (2.22) yields 

Gns = ^~N . (3.2) 

The conductance of the NS junction is quantized in units of Ae 2 /h. This is twice the 
conductance quantum in the normal state, due to the current-doubling effect of Andreev 
reflection [18]. 

In the classical limit N — > oo we recover the well-known result Gns — 2Gn for a classical 
ballistic point contact [12,16,19]. In the quantum regime, however, the simple factor-of-two 
enhancement only holds for the conductance plateaus, where eq. (3.1) applies, and not to the 
transition region between two subsequent plateaus of quantized conductance. To illustrate 
this, we compare in fig. 3 the conductances Gns and 2Gn for Buttiker's model [20] of a 
saddle-point constriction in a two-dimensional electron gas. Appreciable differences appear 
in the transition region, where Gns lies below twice Gn- This is actually a rigorous inequality, 
which follows from eqs. (2.22) and (2.23) for arbitrary transmission matrix: 

G NS < 2G N , V*. (3.3) 
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B. Quantum dot 



Consider next a small confined region (of dimensions comparable to the Fermi wave- 
length), which is weakly coupled by tunnel barriers to two electron reservoirs. We assume 
that transport through this quantum dot occurs via resonant tunneling through a single 
bound state. Let 5 res be the energy of the resonant level, relative to the Fermi level in the 
reservoirs, and let 71 jfi and 72/fr be the tunnel rates through the two barriers. We denote 
7 = 7x + 72 . If 7 <§; AE (with AE the level spacing in the quantum dot), the conductance 
G N in the case of non-interacting electrons has the form 

9^ Gn = F 2 7 + L2 = Tb W» ( 3 - 4 ) 

with Tbw the Breit-Wigner transmission probability at the Fermi level. The normal-state 
transmission matrix ti 2 (e) which yields this conductance has matrix elements [21] 

t 12 ( £ ) = U x t{e)U 2 , r(e) nm = v^mT^ , (3.5) 

£ ^-rcs T 2'7 

where J2 n lin = 7i, J2 n j2n = 72, and U±, U 2 are two unitary matrices (which need not be 
further specified). 

Let us now investigate how the conductance (3.4) is modified if one of the two reservoirs 
is in the superconducting state. The transmission matrix product t 12 t\ 2 (evaluated at the 
Fermi level e — 0) following from eq. (3.5) is 

t 12 t\ 2 = U.MUl M nm = ^ v / 7 t^;. (3.6) 

71 



Its eigenvalues are 



/ ■ 1 ' l>\\ ii ti 1, /,m 
±n ~ l n if 2 <n<N. {6J> 




Substitution into eq. (2.22) yields the conductance 

r 4e ' ( Tbw 2 - 4e ' ( 27172 ^ 2 f 3 R) 

GNS ~ ~h {2-I- ) ~ X U 2 cs + 7 2 + 7 2 2 J • (3 ' 8) 

The conductance on resonance (e res = 0) is maximal in the case of equal tunnel rates 
(71 — 72), and is then equal to 4e 2 /h — independent of 7. The lineshape for this case is 
shown in fig. 4 (solid curve). It differs substantially from the Lorentzian lineshape (3.4) of 
the Breit-Wigner formula (dotted curve). 

The amplitude and lineshape of the conductance resonance (3.8) does not depend on 
the relative magnitude of the resonance width 7 and the superconducting energy gap A . 
This is in contrast to the supercurrent resonance in a superconductor - - quantum dot 
- superconductor Josephson junction, which depends sensitively on the ratio 7//A [22,23]. 
The difference can be traced to the fact that the conductance (in the zero-temperature, zero- 
voltage limit) is strictly a Fermi-level property, whereas all states within A of the Fermi 
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FIG. 4. Conductance versus energy of the resonant level, from eq. (3.8) for the case of equal 
tunnel barriers (solid curve). The dotted curve is the Breit-Wigner transmission probability (3.4). 
The inset shows schematically the normal-metal — quantum-dot — superconductor junction. 



level contribute to the Josephson effect. (For an extension of eq. (3.8) to finite voltages, 
see ref. [24].) Since we have assumed non-interacting quasiparticles, the above results apply 
to a quantum dot with a small charging energy U for double occupancy of the resonant 
state. Devyatov and Kupriyanov [25], and Hekking et al. [26], have studied the influence 
of Coulomb repulsion on resonant tunneling through an NS junction, in the temperature 
regime k-&T ^> 7 where the resonance is thermally broadened. The extension to the low- 
temperature regime of an intrinsically broadened resonance remains to be investigated. 



C. Disordered junction 

We now turn to the regime of diffusive transport through a disordered point contact or 
microbridge between a normal and a superconducting reservoir. The model considered is 
that of an NS junction containing a disordered normal region of length L much greater than 
the mean free path I for elastic impurity scattering, but much smaller than the localization 
length Nl. We calculate the average conductance of the junction, averaged over an ensemble 
of impurity configurations. We begin by parameterizing the transmission eigenvalue T n in 
terms of a channel-dependent localization length ( n : 

Tn = cosh 2 (L/C)' (3 - 9) 

A fundamental result in quantum transport is that the inverse localization length is uniformly 
distributed between and 1/Cmin — 1 / 1 for I <^ L <^ Nl [27-30]. One can therefore write 

00 

dxf(cosh- 2 x), (3.10) 



En=l f(Tn)) J L/ ^dxf(cOsh- 2 x) 



En=lTn) J^dx COsh" 2 X 
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where (...) indicates the ensemble average and f(T) is an arbitrary function of the trans- 
mission eigenvalue such that f(T) — > for T — > 0. In the second equality in eq. (3.10) we 
have used that L/( min ~ L/l >> 1 to replace the upper integration limit by oo. 
Combining eqs. (2.22), (2.23), and (3.10), we find 

roo ( posh -2 t \ ^ 
(Gns) =2(G N > I dx [ 2 _ cQsh -2 x j = (Gn>. (3.H) 

We conclude that — although Gns according to eq. (2.22) is of second order in the trans- 
mission eigenvalues T n — the ensemble average (Gns) is of first order ml/L. The resolution 
of this paradox is that the T's are not distributed uniformly, but are either exponentially 
small (closed channels) or of order unity (open channels) [28] . Hence the average of T 2 is of 
the same order as the average of T n . Off-diagonal elements of the transmission matrix ttf 
are crucial to arrive at the result (3.11). Indeed, if one would evaluate eq. (2.22) with the 
transmission eigenvalues T n replaced by the modal transmission probabilities T n , one would 
find a totally wrong result: Since 1J, ~ l/L < 1, one would find Gns — (l/L)Gw — which 
underestimates the conductance of the NS junction by the factor L/l. 

Previous work [31,32] had obtained the equality of Gns and Gn from semiclassical equa- 
tions of motion, as was appropriate for macroscopic systems which are large compared to the 
normal-metal phase- coherence length 1$. The present derivation, in contrast, is fully quan- 
tum mechanical. It applies to the "mesoscopic" regime L < 1$, in which transport is phase 
coherent. Takane and Ebisawa [33] have studied the conductance of a disordered phase- 
coherent NS junction by numerical simulation of a two-dimensional tight-binding model. 
They found (Gns) = (Gn) within numerical accuracy for I <C L <C Nl, in agreement with 
eq. (3.11). 

If the condition L <C Nl is relaxed, differences between (Gns) an d (Gn) appear. To 
lowest order in L/Nl, the difference is a manifestation of the weak-localization effect, as we 
discuss in the following section. 



IV. WEAK LOCALIZATION 

An NS junction shows an enhanced weak-localization effect, in comparison with the 
normal state [15]. The origin of the enhancement can be understood in a simple way, as 
follows. 

We return to the parameterization T n = 1/ cosh 2 (L/£ n ) introduced in eq. (3.9), and 
define the density of localization lengths p((,L) = (J2n$(( — Cn))i- The subscript L refers 
to the length of the disordered region. Using the identity cosh2x = 2 cosh 2 x — 1, the 
ensemble-average of eq. (2.22) becomes 

(Gns)l = -yJ o dCp(C,£)cosh- 2 (2L/C). (4.1) 
In the same parameterization, one has 

(G n )l = X / q dCp(C,£)cosh- 2 (L/C). (4.2) 
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In the "open-channel approximation" [34], the integrals over ( are restricted to the range 
( > L of localization lengths greater than the length of the conductor. In this range 
the density p((, L) is approximately independent of L. The whole L-dependence of the 
integrands in eqs. (4.1) and (4.2) lies then in the argument of the hyperbolic cosine, so that 

(G N s)l = 2(G N ) 2L . (4.3) 

This derivation formalizes the intuitive notion that Andreev reflection at an NS interface 
effectively doubles the length of the normal- metal conductor [33]. 

Consider now the geometry W <C L relevant for a microbridge. In the normal state one 

has 

(G N > = (W7L)<7 Drude - 5G N , (4.4) 

where 0"D ru de is the classical Drude conductivity. The L-independent term 5G^ is the weak- 
localization correction, given by [35] 5Gn = § e 2 '/h. Equation (4.3) then implies that 

(G N s> = (W/L)a Drude - 5G NS , (4.5) 

with <5Gns = 2 5G N . We conclude that Andreev reflection increases the weak- localization 
correction, by a factor of two according to this qualitative argument [15]. A rigorous the- 
ory [36-38] of weak localization in an NS microbridge shows that the increase is actually 
somewhat less than a factor of two, 1 

SG m = {2-8n- 2 )e 2 /h= 1.78 5G N . (4.6) 

As pointed out in ref. [39], the enhancement of weak localization in an NS junction can be 
observed experimentally as a dip in the differential conductance (jns(^0 — 01/ dV around 
zero voltage. The dip occurs because an applied voltage destroys the enhancement of weak 
localization by Andreev reflection, thereby increasing the conductance by an amount 

5G m - SGx « 0.5 e 2 /h (4.7) 

at zero temperature. [At finite temperatures, we expect a reduction of the size of dip by a 
factor 2 (L c /L) 2 , where L c = min (1$, JfiD/k^T) is the length over which electrons and holes 



1 Equation (4.6) follows from the general formula 5 A = \a(l) + f£°dx (4x 2 + ir 2 ) 1 a(cosh 2 x) for 
the weak-localization correction in a wire geometry, where A is an arbitrary transport property of 
the form A = E„a(^n)- 

2 The reduction factor (L c /L) 2 for the size of the conductance dip when W < L c < L is estimated 
as follows: Consider the wire as consisting of L/L c phase-coherent segments of length L c in series. 
The first segment, adjacent to the superconductor, has a conductance dip 5G\ ~ e 2 /h, while the 
other segments have no conductance dip. The resistance R\ of a single segment is a fraction L c /L 
of the total resistance R of the wire. Since 5G/G = —5R/R = —5R\/R and 5R\ = —R 2 5G\ ~ 
-{L c /L) 2 R 2 e 2 /h, we find SG ~ (L c /L) 2 e 2 /h. 
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remain phase coherent.] We emphasize that in the normal state, weak localization can not 
be detected in the current-voltage characteristic. The reason why a dip occurs in Gns(^) 
and not in G^(V) is that an applied voltage (in contrast to a magnetic field) does not break 
time-reversal symmetry — but only affects the phase coherence between the electrons and 
the Andreev-reflected holes (which differ in energy by up to 2eV). The width V c of the 
conductance dip is of the order of the Thouless energy E c = nhD/L 2 (with D the diffusion 
coefficient of the junction; L should be replaced by L c if L > L c ). This energy scale is such 
that an electron and a hole acquire a phase difference of order tt on traversing the junction. 
The energy E c is much smaller than the superconducting energy gap A , provided L £ 
(with £ ~ (hD/Ao) 1 / 2 the superconducting coherence length in the dirty-metal limit). The 
separation of energy scales is important, in order to be able to distinguish experimentally 
the current due to Andreev reflection below the energy gap from the quasi-particle current 
above the energy gap. 

The first measurement of the conductance dip predicted in ref. [39] has been reported 
recently by Lenssen et al. [40]. The system studied consists of the two-dimensional electron 
gas in a GaAs/AlGaAs heterostructure with Sn/Ti superconducting contacts (W = 10 /xm, 
L = 0.8 jum). No supercurrent is observed, presumably because 1$ ~ 0.4 /im is smaller 
than L. (The phase-coherence length 1$ is estimated from a conventional weak-localization 
measurement in a magnetic field.) The data for the differential conductance is reproduced 
in fig. 5. At the lowest temperatures (10 mK) a rather small and narrow conductance 
dip develops, superimposed on a large and broad conductance minimum. The size of the 
conductance dip is about 2 e 2 /h. Since in the experimental geometry W > L > 1$, and there 
are two NS interfaces, we would expect a dip of order 2(W/l l j > )(l < p/ L) 2 x 0.5e 2 /h ~ 6e 2 /h, 
simply by counting the number of phase-coherent segments adjacent to the superconductor. 
This is three times as large as observed, but the presence of a tunnel barrier at the NS 
interface might easily account for this discrepancy. (The Schottky barrier at the interface 
between a semiconductor and superconductor presents a natural origin for such a barrier.) 
The conductance dip has width V c ~ 0.25 mV, which is less than the energy gap A = 
0.56 meV of bulk Sn — but not by much. Experiments with a larger separation of energy 
scales are required for a completely unambiguous identification of the phenomenon. 

An essential requirement for the appearance of a dip in the differential conductance is 
a high probability for Andreev reflection at the NS boundary. This is illustrated in fig. 
6, which shows the results of numerical simulations [39] of transport through a disordered 
normal region connected via a tunnel barrier to a superconductor. The tunnel barrier is 
characterized by a transmission probability per mode T. The dash-dotted lines refer to an 
ideal interface (r = 1), and show the conductance dip due to weak localization, discussed 
above. For T ~ 0.2-0.4 the data for Gns (filled circles) shows a crossover 3 to a conductance 
peak. This is the phenomenon of reflectionless tunneling, discussed in the following section. 



3 The crossover is accompanied by an "overshoot" around eV & E c , indicating the absence of an 
"excess current" (i.e. the linear I-V characteristic for eV S> E c extrapolates back through the 
origin). We do not have an analytical explanation for the overshoot. 
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FIG. 5. Differential conductance as a function of applied voltage at three different temperatures. 
Experimental data by Lenssen et al. for a two-dimensional electron gas with superconducting 
contacts. The dip around zero voltage, which is superimposed on the broad minimum at the lowest 
temperature, is attributed to the enhancement of weak localization by Andreev reflection. (From 
ref. [40].) 

V. REFLECTIONLESS TUNNELING 

In 1991, Kastalsky et al. [41] discovered a large and narrow peak in the differential 
conductance of a Nb-InGaAs junction. We reproduce their data in fig. 7. (A similar peak 
is observed as a function of magnetic field.) Since then a great deal of experimental [42-48], 
numerical [39,49], and analytical work [50-54] has been done on this effect. Here we focus 
on the explanation in terms of disorder-induced opening of tunneling channels [30,54], which 
is the most natural from the view point of the scattering formula (2.22), and which we feel 
captures the essence of the effect. Equivalently, the conductance peak can be explained in 
terms of a non-equilibrium proximity effect, which is the preferred explanation in a Green's 
function formulation of the problem [52,55-57]. We begin by reviewing the numerical work 
[39]. 

A. Numerical simulations 

A sharp peak in the conductance around V, B = is evident in the numerical simulations 
for T = 0.2 (dotted lines in fig. 6). While Gn depends only weakly on B and V in this range 
(open circles), Gns drops abruptly (filled circles). The width of the conductance peak in B 
and eV is respectively of order B c = h/eLW (one flux quantum through the normal region) 
and eV c = irhD/L 2 = E c (the Thouless energy). The width of the peak is the same as 
the width of the conductance dip due to weak localization, which occurs for larger barrier 
transparencies. The size of the peak is much greater than the dip, however. 

It is instructive to first discuss the classical resistance -Rns ss °f the NS junction. The 
basic approximation in R^ s is that currents rather than amplitudes are matched at the NS 
interface [31]. The result is 
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FIG. 6. Voltage and magnetic field dependence of Gns (filled circles) and Gn (open circles). 
Numerical simulation of a disordered normal region (L/W = 4.8, l/L = 0.12, N = 14) in series 
with a tunnel barrier (transmission probability per mode T; dotted lines: T = 0.2; dashed: T = 0.6; 
dash-dotted: T = 1). Note the crossover from a dip (weak localization) to a peak (reflectionless 
tunneling) in Gns on reducing T. (From ref. [39].) 



= (h/2Ne 2 ) [L/l + 2T' 2 + 0(1)\ . (5.1) 

The contribution from the barrier is oc T~ 2 because tunneling into a superconductor is a 
two-particle process [58]: Both the incident electron and the Andreev-reflected hole have to 
tunnel through the barrier (the net result being the addition of a Cooper pair to the super- 
conducting condensate [1]). Equation (5.1) is to be contrasted with the classical resistance 
^>ciass j n norma } state, 

R^ s = (h/2Ne 2 ) [L/l + I^ 1 + 0(1)] , (5.2) 

where the contribution of a resistive barrier is oc In the absence of a tunnel barrier (i.e. 
for T = 1), -Rns ss = -^N ass f° r L ^> I, in agreement with refs. [31,32]. Let us now see how 
these classical results compare with the simulations [39] . 

In fig. 8 we show the resistance (at V = 0) as a function of T in the absence and presence of 
a magnetic field. (The parameters of the disordered region are the same as for fig. 6.) There 
is good agreement with the classical eqs. (5.1) and (5.2) for a magnetic field corresponding to 
10 flux quanta through the disordered segment (fig. 8b). For B = 0, however, the situation 
is different (fig. 8a). The normal-state resistance (open circles) still follows approximately 
the classical formula (solid curve). (Deviations due to weak localization are noticeable, but 
small on the scale of the figure.) In contrast, the resistance of the NS junction (filled circles) 
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FIG. 7. Differential conductance (normalized by the normal-state resistance = 0.27 fi) as 
a function of applied voltage at seven different temperatures. Experimental data by Kastalsky et 
al. for a Nb-InGaAs junction. Note the difference with fig. 5: A peak rather than a dip develops 
at the lowest temperatures, and the size of the peak (0.6 f2 _1 ~ 1.5 • 10 4 e 2 //i) is four orders of 
magnitude greater. The width of the peak is comparable to the width of the dip in fig. 5. (From 
ref. [41].) 



lies much below the classical prediction (dotted curve). The numerical data shows that for 
r ^> l/L one has approximately 

R NS (B = 0,V = 0)aB$r, (5.3) 

which for T 1 is much smaller than -R^fT 3 - This is the phenomenon of reflectionless 
tunneling: In fig. 8a the barrier contributes to /?ns 111 order just as for single-particle 
tunneling, and not in order T~ 2 , as expected for two-particle tunneling. It is as if the 
Andreev-reflected hole is not reflected by the barrier. The interfering trajectories responsible 
for this effect were first identified by Van Wees et al. [50]. The numerical data of fig. 8a is 
in good agreement with the Green's function calculation of Volkov, Zaitsev, and Klapwijk 
[52] (dashed curve). Both these papers have played a crucial role in the understanding of 
the effect. The scaling theory reviewed below [54] is essentially equivalent to the Green's 
function calculation, but has the advantage of explicitly demonstrating how the opening of 
tunneling channels on increasing the length L of the disordered region induces a transition 
from a T~ 2 dependence to a dependence when L ~ l/Y. 



B. Scaling theory 

We use the parameterization 

cosh 2 x n ' ^ 

similar to eq. (3.9), but now with a dimensionless variable x n G [0, oo). The density of the 
re-variables, for a length L of disordered region, is denoted by 
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FIG. 8. Filled circles: Numerically calculated resistance i?NS of a disordered NS junction, versus 
the transmission probability per mode T of the tunnel barrier at the NS interface; Open circles: 
Resistance i?N of the same junction in the normal state; (a) is for zero magnetic field, (b) is for 
a flux of 10/i/e through the disordered region. The dotted and solid curves are the classical eqs. 
(5.1) and (5.2). The dashed curve is the theory of ref. [52], which for T ^ l/L ~ 0.12 coincides 
with eq. (5.3). (From ref. [39].) 



p(x, L) = (Ej{x ~ Xu))l- (5.5) 

For L = 0, i.e. in the absence of disorder, we have the initial condition imposed by the 
barrier, 

p{x,0) = N5(x-x ), (5.6) 

with T = 1/ cosh 2 xo. The scaling theory describes how p(x,L) evolves with increasing L. 
This evolution is governed by the equation 



\ f°° 

— p(x, s) = - — — p(x,s)— dx' p(x , ,s)\n\sinh 2 x-smh 2 x'\, (5.7) 



where we have defined s = L/l. This non-linear diffusion equation was derived by Mello and 
Pichard [59] from a Fokker-Planck equation [34,60,61] for the joint distribution function of 
all N eigenvalues, by integrating out N — 1 eigenvalues and taking the large- N limit. This 
limit restricts its validity to the metallic regime (N ^> L/l), and is sufficient to determine 
the leading order contribution to the average conductance, which is O(N). The weak- 
localization correction, which is 0(1), is neglected here. A priori, eq. (5.7) holds only for a 
"quasi-one-dimensional" wire geometry (length L much greater than width W) , because the 
Fokker-Planck equation from which it is derived requires L ^> W. Numerical simulations 
indicate that the geometry dependence only appears in the 0(1) corrections, and that the 
O(N) contributions are essentially the same for a wire, square, or cube. 

In ref. [54] it is shown how the scaling equation (5.7) can be solved exactly, for arbitrary 
initial condition p(x, 0) = po(x). The method of solution is based on a mapping of eq. (5.7) 
onto Euler's equation for the isobaric flow of a two-dimensional ideal fluid: L corresponds 
to time and p to the y-component of the velocity field on the x-axis. [Please note that in 
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FIG. 9. Eigenvalue density p(x, s) as a function of x (in units of s = L/l) for V = 0.1. Curves 
a,b,c,d,e are for s = 2,4,9,30,100, respectively. The solid curves are from eq. (5.8), the dashed 
curves from eq. (5.12). The collision of the density profile with the boundary at x = 0, for 
s = s c = (1 — r)/r, signals the disorder-induced opening of tunneling channels responsible for the 
reflectionless-tunneling effect. (From ref. [54]-) 



this section x is the auxiliary variable defined in eq. (5.4) and not the physical coordinate 
in fig. 2.] The result is 



p(x, s) = (2N/n) ImC/(x - i0 + , s), 
where the complex function U (z, s) is determined by 

U(z,s) = U (z-sU(z,s)). 
The function U (z) is fixed by the initial condition, 

sinh 2z f°° , , Po(x') 



U (z) = 



2N 



dx' 



sinh z — sinh x' 



(5.8) 



(5.9) 



(5.10) 



The implicit equation (5.9) has multiple solutions in the entire complex plane; We need 
the solution for which both z and z — sU (z, s) lie in the strip between the lines y = and 
y = — 7r/2, where z = x + iy. 

The initial condition (5.6) corresponds to 



U (z) = i sinh 2,2 (cosh 2 2 -r^ 1 )- 1 . 



(5.11) 



The resulting density (5.8) is plotted in fig. 9 (solid curves), for T = 0.1 and several values 
of s. For s ^> 1 and x <^ s it simplifies to 



x = |arccoshr — |Ts(r 2 — l) 1 / 2 cosa, 
a = irsN^ 1 p(x, s), r = a(Ts sin a) -1 , 



(5.12) 



shown dashed in fig. 9. Equation (5.12) agrees with the result of a Green's function cal- 
culation by Nazarov [30]. For s = (no disorder), p is a delta function at x . On adding 
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disorder the eigenvalue density rapidly spreads along the x-axis (curve a), such that p < N/ s 
for s > 0. The sharp edges of the density profile, so uncharacteristic for a diffusion profile, 
reveal the hydrodynamic nature of the scaling equation (5.7). The upper edge is at 

x ma x = s + |ln( S /r) + 0(l). (5.13) 

Since L/x has the physical significance of a localization length [34], this upper edge corre- 
sponds to a minimum localization length £ min = L/x max of order /. The lower edge at x niin 
propagates from x to in a "time" s c — (1 — T)/T. For 1 < s < s c one has 

x min = iarccosh (s c /s) - \[1 - (s/s,) 2 } 1 ' 2 . (5.14) 

It follows that the maximum localization length £ max = L/:r mm increases if disorder is 
added to a tunnel junction. This paradoxical result, that disorder enhances transmission, 
becomes intuitively obvious from the hydrodynamic correspondence, which implies that 
p(x, s) spreads both to larger and smaller x as the fictitious time s progresses. When s = s c 
the diffusion profile hits the boundary at x = (curve c), so that x min = 0. This implies 
that for s > s c there exist scattering states (eigenfunctions of ttf) which tunnel through 
the barrier with near-unit transmission probability, even if T <C 1. The number N opcn of 
transmission eigenvalues close to one (open channels) is of the order of the number of x n 's 
in the range to 1 (since T n = 1/ cosh 2 x„ vanishes exponentially if x n > 1). For s ^> s c 
(curve e) we estimate 

N opcn ~ p(0, s) = N(s + r- 1 )" 1 , (5.15) 

where we have used eq. (5.12). The disorder-induced opening of tunneling channels was 
discovered by Nazarov [30]. It is the fundamental mechanism for the T~ 2 to T _1 transition 
in the conductance of an NS junction, as we now discuss. 

According to eqs. (2.22), (2.23), (5.4), and (5.5), the average conductances (Gns) an d 
(Gn) are given by the integrals 

(Gns) = —r- I drr p(x, s) cosh -2 2x, (5.16) 
h Jo 

(Cn) = —r~ / dx p(x, s) cosh~ 2 x. (5-17) 
h Jo 

Here we have used the same trigonometric identity as in eq. (4.1). For T ^> l/L one is in 
the regime s ^> s c of curve e in fig. 9. Then the dominant contribution to the integrals 
comes from the range x/s <C 1 where p(x,s) ~ p(0, s) = N(s + T -1 ) -1 is approximately 
independent of x. Substitution of p(x,s) by p(0, s) in eqs. (5.16) and (5.17) yields directly 

(Gns) ~ (C N ) ~ l/^ ass , (5.18) 

in agreement with the result (5.3) of the numerical simulations. 

Equation (5.18) has the linear T dependence characteristic for reflectionless tunneling. 
The crossover to the quadratic T dependence when r < l/L is obtained by evaluating the 
integrals (5.16) and (5.17) with the density p(x, s) given by eq. (5.8). The result is [54] 
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FIG. 10. Dependence of the resistance i?NS on the length L of the disordered normal region 
(hatched in the inset), for different values of the transmittance V of the NS interface. Solid curves 
are computed from eq. (5.19), for T = 1, 0.8, 0.6, 0.4, 0.1 from bottom to top. For r<l the dashed 
curve is approached. (From ref. [54]-) 



(G m ) = (2Ne 2 /h)(s + Q- 1 )-\ 
(G N ) = (2Ne 2 /h)(s + T- 1 )- 1 - 

The "effective" tunnel probability Q is denned by 

e ( e 



Q 



s cos 9 \ Ts cos 9 



1 + sin#) - 1 



where 9 G (0, ir/2) is the solution of the transcendental equation 

9[1 - ±r(l -sin0)] = rscosfl. 



(5.19) 
(5.20) 



(5.21) 



(5.22) 



For T < 1 (or s > 1) eqs. (5.21) and (5.22) simplify to Q = Tsin9, 9 = Tscos9, in 
precise agreement with the Green's function calculation of Volkov, Zaitsev, and Klapwijk 
[52]. According to eq. (5.20), the normal-state resistance increases linearly with the length 
L of the disordered region, as expected from Ohm's law. This classical reasoning fails if one 
of the contacts is in the superconducting state. The scaling of the resistance i?NS = l/(G r Ns) 
with length, computed from eq. (5.19), is plotted in fig. 10. For T = 1 the resistance increases 
monotonically with L. The ballistic limit L — > equals h/4Ne 2 , half the contact resistance 
of a normal junction because of Andreev reflection (cf. section 3.1). For T < 0.5 a resistance 
minimum develops, somewhat below L = l/Y. The resistance minimum is associated with 
the crossover from a quadratic to a linear dependence of i?NS on 1/r. 

If Ts 3> 1 one has 9 — > n/2, hence Q — > T. In the opposite regime Ts <^ 1 one has 9 — > Ts, 
hence Q — ^ T 2 s. The corresponding asymptotic expressions for (G NS ) are (assuming T 1 
and s> 1): 



(Gns) = (2Ne 2 /h)(s + T~ l )-\ if Ts > 1, 
(Gns) = (2Ne 2 /h)T 2 s, if Ts < 1. 



(5.23) 
(5.24) 
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In either limit the conductance is greater than the classical result 



G^ ss = (2Ne 2 /h)(s + 2T- 2 )-\ (5.25) 

which holds if phase coherence between electrons and holes is destroyed by a voltage or 
magnetic field. The peak in the conductance around V, B = is of order AG NS = (Gns) — 
Gn| ss , which has the relative magnitude 

<GW " 2 + P? (5 ' 26) 

The scaling theory assumes zero temperature. Hekking and Nazarov [53] have studied 
the conductance of a resistive NS interface at finite temperatures, when L is greater than 



the correlation length L c = min (l^, \JhD /ksT). Their result is consistent with the limiting 
expression (5.24), if s = L/l is replaced by L c /l. The implication is that, if L > L c , the 
non-linear scaling of the resistance shown in fig. 10 only applies to a disordered segment 
of length L c adjacent to the superconductor. For the total resistance one should add the 
Ohmic contribution of order (h/e 2 )(L — L c )/l from the rest of the wire. 



C. Double-barrier junction 

In the previous subsection we have discussed how the opening of tunneling channels (i.e. 
the appearance of transmission eigenvalues close to one) by disorder leads to a minimum 
in the resistance when L ~ l/T. The minimum separates a from a T~ 2 dependence 
of the resistance on the transparency of the interface. We referred to the T" 1 dependence 
as "reflectionless tunneling", since it is as if one of the two quasiparticles which form the 
Cooper pair can tunnel through the barrier with probability one. In the present subsection 
we will show, following ref . [62] , that a qualitatively similar effect occurs if the disorder in the 
normal region is replaced by a second tunnel barrier (tunnel probability T') . The resistance 
at fixed r shows a minimum as a function of I~" when T' ~ T. For I~" < T the resistance has 
a r _1 dependence, so that we can speak again of reflectionless tunneling. 

We consider an NIiNI 2 S junction, where N = normal metal, S = superconductor, and 
Ij = insulator or tunnel barrier (transmission probability per mode 1^ = 1/ cosh 2 a:i). We 
assume ballistic motion between the barriers. (The effect of disorder is discussed later.) A 
straightforward calculation yields the transmission probabilities T n of the two barriers in 
series, 

T n = (a + bcoscpn)- 1 , (5.27) 
a = \ + \ cosh2«i cosh2o; 2 , b = \ sinh 2a\ sinh 2a 2 , (5.28) 

where f n is the phase accumulated between the barriers by mode n. Since the transmission 
matrix t is diagonal, the transmission probabilities T n are identical to the eigenvalues of tt^. 
We assume that L ^> Ap (Ap is the Fermi wavelength) and iVTj ^> 1, so that the conductance 
is not dominated by a single resonance. In this case, the phases tp n are distributed uniformly 
in the interval (0, 27r) and we may replace the sum over the transmission eigenvalues in eqs. 
(2.22) and (2.23) by integrals over <p: E„=i f(<p n ) -> (N/2rt) J 27T dip f(tp). The result is 
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FIG. 11. Dependence of the resistances i?N and i?NS of ballistic NININ and NINIS structures, 
respectively, on barrier transparency T±, while transparency r 2 = 0.1 is kept fixed [computed from 
eqs. (5.29) and (5.30)]. The inset shows the NINIS structure considered. (From ref. [62].) 
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(5.29) 
(5.30) 



These expressions are symmetric in the indices 1 and 2: It does not matter which of the 
two barriers is closest to the superconductor. In the same way we can compute the entire 
distribution of the transmission eigenvalues, p(T) = I] n 5(T — T n ) — > (N/2tt) J 27r dcpS(T — 
T((p)). Substituting T(ip) = (a + bcosip)^ 1 from eq. (5.27), one finds 



p(T) = ^ (6 2 T 2 - (aT - I) 2 ) 



-1/2 



TFT 

In fig. 1 1 we plot the resistance = 1 / and R^s 
and (5.30). Notice that R^ follows Ohm's law, 



(5.31) 



1/GnS) following from eqs. (5.29) 



Rn — 



h 



2Ne 2 



(i/r 1 + i/r 2 -i), 



(5.32) 



as expected from classical considerations. In contrast, the resistance i?NS has a minimum 
if one of the T's is varied while keeping the other fixed. This resistance minimum cannot 
be explained by classical series addition of barrier resistances. If T 2 1 is fixed and I\ is 
varied, as in fig. 11, the minimum occurs when I\ = v^2T 2 . The minimal resistance R^s 1 is 
of the same order of magnitude as the resistance R^ in the normal state at the same value 
of Ti and T 2 . In particular, we find that -R^g 1 depends linearly on 1/Tj, whereas for a single 
barrier R NS oc 1/T 2 . 
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FIG. 12. Density of transmission eigenvalues through a normal region containing a potential 
barrier (transmission probability T = 0.4). The left panel (a) shows the disorder-induced opening 
of tunneling channels (solid curve: s = 0.04; dotted: s = 0.4; dashed: s = 5; where s = L/l). The 
right panel (b) shows the opening of channels by a second tunnel barrier (transparency T'; solid 
curve: V = 0.95; dotted: T' = 0.8; dashed: V = 0.4). The curves in (a) are computed from eq. 
(5.8), the curves in (b) from eq. (5.31). (From ref. [62].) 



The linear dependence on the barrier transparency shows the qualitative similarity of a 
ballistic NINIS junction to the disordered NIS junction considered in the previous subsection. 
To illustrate the similarity, we compare in fig. 12 the densities of normal-state transmission 
eigenvalues. The left panel is for an NIS junction [computed using eq. (5.8)], the right panel 
is for an NINIS junction [computed from eq. (5.31)]. In the NIS junction, disorder leads to 
a bimodal distribution p(T), with a peak near zero transmission and another peak near unit 
transmisssion (dashed curve). A similar bimodal distribution appears in the ballistic NINIS 
junction, for approximately equal transmission probabilities of the two barriers. There are 
also differences between the two cases: The NIS junction has a uni-modal p(T) if L/l < 1/T, 
while the NINIS junction has a bimodal p(T) for any ratio of Y\ and r 2 . In both cases, 
the opening of tunneling channels, i.e. the appearance of a peak in p(T) near T = 1, is the 
origin for the 1 /Y dependence of the resistance. 

The scaling equation of section 5.2 can be used to investigate what happens to the 
resistance minimum if the region of length L between the tunnel barriers contains impurities, 
with elastic mean free path I. As shown in ref. [62], the resistance minimum persists as long 
as I ^ TL. In the diffusive regime (I <C L) the scaling theory is found to agree with 
the Green's function calculation by Volkov, Zaitsev, and Klapwijk for a disordered NINIS 
junction [52]. For strong barriers (Ti, r 2 C 1) and strong disorder (L I), one has the two 
asymptotic formulas 

2Ne 2 rjrl 



Gns — 



Gns — 



h (r? + r| 

2Ne 2 



.3/2 : 



if ri,r 2 < i/l, 



h 



(L/i + i/r 1 + i/r 2 )-\ if r 1 ,r 2 »//L. 



(5.33) 



(5.34) 



Equation (5.33) coincides with eq. (5.29) in the limit a± : a 2 ^> 1 (recall that 1^ = 1/ cosh 2 otj). 
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This shows that the effect of disorder on the resistance minimum can be neglected as long 
as the resistance of the junction is dominated by the barriers. In this case Gns depends 
linearly on F 1 and T 2 only if F 1 pa T 2 - Equation (5.34) shows that if the disorder dominates, 
Gns has a linear T-dependence regardless of the relative magnitude of I\ and T 2 . 

We have assumed zero temperature, zero magnetic field, and infinitesimal applied voltage. 
Each of these quantities is capable of destroying the phase coherence between the electrons 
and the Andreev-reflected holes, which is responsible for the resistance minimum. As far 
as the temperature T and voltage V are concerned, we require ksT, eV <C h/r dwe n for the 
appearance of a resistance minimum, where Td we ii is the dwell time of an electron in the 
region between the two barriers. For a ballistic NINIS junction r dwcll ~ L/v F T, while for a 
disordered junction Td wc ii — L 2 /vfTI is larger by a factor L/l. It follows that the condition 
on temperature and voltage becomes more restrictive if the disorder increases, even if the 
resistance remains dominated by the barriers. As far as the magnetic field B is concerned, 
we require B <C h/eS (with S the area of the junction perpendicular to B), if the motion 
between the barriers is diffusive. For ballistic motion the trajectories enclose no flux, so no 
magnetic field dependence is expected. 

A possible experiment to verify these results might be scanning tunneling microscopy 
(STM) of a metal particle on a superconducting substrate [63]. The metal-superconductor 
interface has a fixed tunnel probability T 2 . The probability Ti for an electron to tunnel from 
STM to particle can be controlled by varying the distance. (Volkov has recently analyzed 
this geometry in the regime that the motion from STM to particle is diffusive rather than by 
tunneling [64].) Another possibility is to create an NINIS junction using a two-dimensional 
electron gas in contact with a superconductor. An adjustable tunnel barrier could then be 
implemented by means of a gate electrode. 

D. Circuit theory 

The scaling theory of ref. [54], which was the subject of section 5.2, describes the tran- 
sition from the ballistic to the diffusive regime. In the diffusive regime it is equivalent to 
the Green's function theory of ref. [52]. A third, equivalent, theory for the diffusive regime 
was presented recently by Nazarov [65]. Starting from a continuity equation for the Keldysh 
Green's function [66], and applying the appropriate boundary conditions [67], Nazarov was 
able to formulate a set of rules which reduce the problem of computing the resistance of 
an NS junction to a simple exercise in circuit theory. Furthermore, the approach can be 
applied without further complications to multi-terminal networks involving several normal 
and superconducting reservoirs. Because of its practical importance, we discuss Nazarov's 
circuit theory in some detail. 

The superconductors Si should all be at the same voltage, but may have a different 
phase (pi of the pair potential. Zero temperature is assumed, as well as infinitesimal voltage 
differences between the normal reservoirs (linear response). The reservoirs are connected by 
a set of diffusive normal-state conductors (length L iy mean free path If = Li/U ^> 1). 
Between the conductors there may be tunnel barriers (tunnel probability Ti). The presence 
of superconducting reservoirs has no effect on the resistance (h/2Ne 2 )si of the diffusive 
conductors, but affects only the resistance h/2Ne 2 Vf of the tunnel barriers. The tunnel 



26 



probability Tj of barrier i is renormalized to an effective tunnel probability rf ff , which 
depends on the entire circuit. 

Nazarov's rules to compute the effective tunnel probabilities are as follows. To each 
node and to each terminal of the circuit one assigns a vector Hi of unit length. For a 
normal reservoir, Hi = (0, 0, 1) is at the north pole, for a superconducting reservoir, Hi = 
(cos0j, sin0j, 0) is at the equator. For a node, Hi is somewhere on the northern hemisphere. 
The vector Hi is called a "spectral vector" , because it is a particular parameterization of the 
local energy spectrum. If the tunnel barrier is located between spectral vectors Hi and n 2 , 
its effective tunnel probability is 4 

T cfr = (ni • n 2 )r = T cos 9 12 , (5.35) 

where 9 12 is the angle between Hi and n 2 . The rule to compute the spectral vector of node 
% follows from the continuity equation for the Green's function. Let the index k label the 
nodes or terminals connected to node % by a single tunnel barrier (with tunnel probability 
r&). Let the index q label the nodes or terminals connected to i by a diffusive conductor 
(with L/l = s q ). The spectral vectors then satisfy the sum rule [65] 

x H k )T k + ]T(n, x H q ) y 1 q> = 0. (5.36) 

k q Sq^Jl - (Hi ■ Hg) 2 

This is a sum rule for a set of vectors perpendicular to Hi of magnitude T^sin^ or 9i q /s q , 
depending on whether the element connected to node % is a tunnel barrier or a diffusive 
conductor. There is a sum rule for each node, and together the sum rules determine the 
spectral vectors of the nodes. 

As a simple example, let us consider the system of section 5.2, consisting of one normal 
terminal (N), one superconducting terminal (S), one node (labeled A), and two elements: A 
diffusive conductor (with L/l = s) between N and A, and a tunnel barrier (tunnel probability 
T) between A and S (see fig. 13). There are three spectral vectors, %, H s , and n\. All 
spectral vectors lie in one plane. (This holds for any network with a single superconducting 
terminal.) The resistance of the circuit is given by R = (h/2Ne 2 )(s + 1/T cff ), with the 
effective tunnel probability 

r eff = Tcos^as = rsintf. (5.37) 

Here 9 e [0, n/2] is the polar angle of n\. This angle is determined by the sum rule (5.36), 
which in this case takes the form 

Tcos9 -9/s = 0. (5.38) 

Comparison with section 5.2 shows that r eff coincides with the effective tunnel probability 
Q of eq. (5.21) in the limit s> 1, i.e. if one restricts oneself to the diffusive regime. That 
is the basic requirement for the application of the circuit theory. 



4 It may happen that cos 612 < 0, in which case the effective tunnel probability is negative. Nazarov 
has given an example of a four-terminal circuit with T eS < 0, so that the current through this barrier 
flows in the direction opposite to the voltage drop [68]. 
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FIG. 13. At left: Circuit containing two terminals (open circles), one node (filled circle), and 
two elements: A diffusive conductor (shaded) and a tunnel barrier (black). At right: Spectral 
vectors associated with the terminals N,S and with the node A. 
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FIG. 14. Circuit diagram and spectral vectors for a structure containing one normal and two 
superconducting terminals (phase difference 0). 



Let us now consider the "fork junction" of fig. 14, with one normal terminal (N) and 
two superconducting terminals Si and S2 (phases 0i = —0/2 and 2 = 0/2). There is one 
node (A), which is connected to N by a diffusive conductor [Ljl = s), and to Si and S2 
by tunnel barriers (Ti and T 2 ). This structure was studied theoretically by Hekking and 
Nazarov [53] and experimentally by Pothier et al. [69]. For simplicity, let us assume two 
identical tunnel barriers Ti = T 2 = T. Then the spectral vector % = (sin#, 0, cos0) of node 
A lies symmetrically between the spectral vectors of terminals Si and S 2 . The sum rule 
(5.36) now takes the form 

2r|cosi0|cos#-#/s = O. (5.39) 

Its solution determines the effective tunnel rate r eff = r|cos|0|sin^ of each of the two 
barriers in parallel, and hence the conductance of the fork junction, 

2A^p 2 

G = [s + Kricos^lsinfl)- 1 ]- 1 . (5.40) 
Two limiting cases of eqs. (5.39) and (5.40) are 
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FIG. 15. Conductance of a fork junction as a function of magnetic field, showing the depen- 
dence on the phase difference <p of the superconductor at two tunnel barriers. The circles are 
measurements by Pothier et al. [69] of the current / through a Cu wire connected to an oxidized 
Al fork (normal-state resistance i?N = 1.56 kfi). The applied voltage V is sufficiently low that 
I/V is close to the linear-response conductance. (The amplitude of the oscillations at V = is 
3.94- 10~ 6 ft' 1 , somewhat larger than in the figure.) The solid curve is a cosine fit to the data. The 
offset of maximum conductance from B = is attributed to a small residual field in the cryostat. 
( Courtesy of H. Pothier.) 



For = (and 2T — > V) these expressions reduce to the results (5.23) and (5.24) for an 
NS junction with a single superconducting reservoir. The limit (5.42) agrees with the finite- 
temperature result of Hekking and Nazarov [53], if s is replaced by L c /l and a series resistance 
is added due to the normal segment which is further than a correlation length from the NS 
interfaces. The possibility of a dependence of the conductance on the superconducting phase 
difference was noted also in other theoretical works, for different geometries [70-75]. 

The ^-dependence of the conductance of a fork junction has recently been observed by 
Pothier et al. [69]. Some of their data is reproduced in fig. 15. The conductance of a Cu wire 
attached to an oxidized Al fork oscillates as a function of the applied magnetic field. The 
period corresponds to a flux increment of h/2e through the area enclosed by the fork and the 
wire, and thus to A0 = 2tt. The experiment is in the regime where the junction resistance 
is dominated by the tunnel barriers, as in eq. (5. 42). 5 The metal-oxide tunnel barriers in 
such structures have typically very small transmission probabilities (r ~ 10~ 5 in ref. [69]), 
so that the regime of eq. (5.41) is not easily accessible. Larger T's can be realized by the 



5 Equation (5.42) provides only a qualitative description of the experiment, mainly because the 
motion in the arms of the fork is diffusive rather than ballistic. This is why the conductance 
minima in fig. 15 do not go to zero. A solution of the diffusion equation in the actual experimental 
geometry is required for a quantitative comparison with the theory [69]. 




(5.41) 
(5.42) 
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FIG. 16. Numerical calculation of the variance of the fluctuations in Gn and Gns> as a function 
of the average Gn (+ for B = 0; x for a flux of Wh/e). Dotted lines are the analytical results 
from eqs. (6.1) and (6.3). Note the absence of a factor-of-two reduction in Var£?NS on applying a 
magnetic field. (From ref. [39].) 



Schottky barrier at a semiconductor — superconductor interface. It would be of interest to 
observe the crossover with increasing V to the non-sinusoidal ^-dependence predicted by eq. 
(5.41), as a further test of the theory. 



VI. UNIVERSAL CONDUCTANCE FLUCTUATIONS 

So far we have considered the average of the conductance over an ensemble of impurity 
potentials. In fig. 16 we show results of numerical simulations [39] for the variance of the 
sample-to-sample fluctuations of the conductance, as a function of the average conductance 
in the normal state. A range of parameters L, W, I, N was used to collect this data, in the 
quasi-one-dimensional, metallic, diffusive regime I < W < L < Nl. An ideal NS interface 
was assumed (T = 1). The results for VarGN are as expected theoretically [34,35] for 
"universal conductance fluctuations" (UCF): 

VarG N = JLrV//*) 2 - (6.1) 

The index (3 equals 1 in the presence and 2 in the absence of time-reversal symmetry. The 1/(3 
dependence of Var Gn implies that the variance of the conductance fluctuations is reduced 
by a factor of two upon application of a magnetic field, as observed in the simulation (see 
the two dotted lines in the lower part of fig. 16). The data for VarG^s at -B = shows 
approximately a four-fold increase over VarGN- For B ^ 0, the simulation shows that 
VarG^Ns is essentially unaffected by a time-reversal-symmetry breaking magnetic field. In 
contrast to the situation in the normal state, the theory for UCF in an NS junction is quite 
different for zero and for non-zero magnetic field, as we now discuss. 

In zero magnetic field, the conductance of the NS junction is given by eq. (2.22), which 
is an expression of the form A = I]„a(T n ). Such a quantity A is called a linear statistic on 
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the transmission eigenvalues. The word "linear" refers to the fact that A does not contain 
products of different T„'s. The function a(T) may well depend non-linearly on T, as it 
does for Gns 5 where a(T) is a rational function of T. The Landauer formula (2.23) for the 
normal-state conductance is also a linear statistic, with a(T) oc T. It is a general theorem 
in random-matrix theory [76] that the variance of a linear statistic has a 1/(3 dependence 
on the symmetry index (3. Moreover, the magnitude of the variance is independent of 
the microscopic properties of the system (sample size, degree of disorder). This is Imry's 
fundamental explanation for UCF [28]. 

For a wire geometry, there exists a formula for the variance of an arbitrary linear statistic 
[37,77,78], 

2(3n 2 Jo Jo \ dT J \ dV J 

where x(T) = arccoshT -1 / 2 . In the normal state, substitution of a(T) = (2e 2 /h)T into eq. 
(6.2) reproduces the result (6.1). In the NS junction, substitution of a{T) = (4e 2 //i)T 2 (2 - 
T)~ 2 yields, for the case (3 — 1 of zero magnetic field, 

32 

VarG NS = — (2 - 907r- 4 )(e 2 //i) 2 = 4.30 VarG N . (6.3) 
15 

A factor of four between VarG N s and VarG N was estimated by Takane and Ebisawa [33], 
by an argument similar to that which we described in section 4 for the weak-localization 
correction. (A diagrammatic calculation by the same authors [79] gave a factor of six, 
presumably because only the dominant diagram was included.) The numerical data in fig. 
16 is within 10 % of the theoretical prediction (6.3) (upper dotted line). Similar numerical 
results for VarG N s in zero magnetic field were obtained in refs. [33,80]. 

We conclude that UCF in zero magnetic field is basically the same phenomenon for Gn 
and Gns, because both quantities are linear statistics for (3=1. If time-reversal symmetry 
(TRS) is broken by a magnetic field, the situation is qualitatively different. For Gn, broken 
TRS does not affect the universality of the fluctuations, but merely reduces the variance by 
a factor of two. No such simple behavior is to be expected for Gns, since it is no longer a 
linear statistic for (3 = 2. That is a crucial distinction between eq. (2.20) for Gns and the 
Landauer formula (2.23) for Gn, which remains a linear statistic regardless of whether TRS 
is broken or not. This expectation [15] of an anomalous /3-dependence of Var Gns was borne 
out by numerical simulations [39], which showed that the conductance fluctuations in an 
NS junction without TRS remain independent of disorder, and of approximately the same 
magnitude as in the presence of TRS (compare + and x data points in the upper part of 
fig. 16). An analytical theory remains to be developed. 

VII. SHOT NOISE 

The conductance, which we studied in the previous sections, is the time- averaged current 
/ divided by the applied voltage V . Time-dependent fluctuations SI(t) in the current give 
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additional information on the transport processes. The zero-frequency noise power P is 
defined by 

(•CO 

P = A dt(5I(t)5I(0)). (7.1) 

At zero temperature, the discreteness of the electron charge is the only source of fluctuations 
in time of the current. These fluctuations are known as "shot noise" , to distinguish them from 
the thermal noise at non-zero temperature. A further distinction between the two is that the 
shot-noise power is proportional to the applied voltage, whereas the thermal noise does not 
vanish at V = 0. Shot noise is therefore an intrinsically non-equilibrium phenomenon. If the 
transmission of an elementary charge e can be regarded as a sequence of uncorrelated events, 
then P = 2e\I\ = Pp isson as in a Poisson process. In this section we discuss, following ref. 
[81], the enhancement of shot noise in an NS junction. The enhancement originates from 
the fact that the current in the superconductor is carried by Cooper pairs in units of 2e. 
However, as we will see, a simple factor-of-two enhancement applies only in certain limiting 
cases. 

In the normal state, the shot-noise power (at zero temperature and infinitesimal applied 
voltage) is given by [82] 

N 

P N = P Trttt(i _ tt t) =P J2 T n (l - T n ), (7.2) 

n=l 

with Pq = 2e\V\ (2e 2 /h). Equation (7.2) is the multi-channel generalization of earlier single- 
channel formulas [83,84]. It is a consequence of the Pauli principle that closed (T n = 0) as 
well as open (T n = 1) scattering channels do not fluctuate and therefore give no contribution 
to the shot noise. In the case of a tunnel barrier, all transmission eigenvalues are small 
(T n <C 1, for all n), so that the quadratic terms in eq. (7.2) can be neglected. Then it follows 
from comparison with eq. (2.23) that Pn = 2e|V^|(jN = 2e\I\ = Pp isson- In contrast, for a 
quantum point contact Pn <C Ppoisson- Since on the plateaus of quantized conductance all the 
T n 's are either or 1, the shot noise is expected to be only observable at the steps between 
the plateaus [84]. For a diffusive conductor of length L much longer than the elastic mean 
free path /, the shot noise Pn = |Pp isson is one-third the Poisson noise, as a consequence of 
noiseless open scattering channels [85,86]. 

The analogue of eq. (7.2) for the shot-noise power of an NS junction is [81] 

N 1 6T 2 n — t i 

P NS = 4P Tr She 4 e (l - She 4 e ) = P £ ( ; { _ T ,; > , (7.3) 

n=l \ Z I n) 

where we have used eq. (2.17) (with e = 0) to relate the scattering matrix s hc for Andreev 
reflection to the transmission eigenvalues T n of the normal region. This requires zero mag- 
netic field. As in the normal state, scattering channels which have T n = or T n = 1 do 
not contribute to the shot noise. However, the way in which partially transmitting channels 
contribute is entirely different from the normal state result (7.2). 

Consider first an NS junction without disorder, but with an arbitrary transmission prob- 
ability T per mode of the interface. In the normal state, eq. (7.2) yields P N = (1 — r)Pp isson, 
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implying full Poisson noise for a high tunnel barrier (r <C 1). For the NS junction we find 
from eq. (7.3) 



i6r 2 (i-r) _ 8(1 - r) 
(2 - r) 4 ~ (2 - ry 



Pns = PoN fn 1 ' = f \ ; P Poisson , (7.4) 



where in the second equality we have used eq. (2.22). This agrees with results obtained 
by Khlus [83], and by Muzykantskii and Khmel'nitskii [87], using different methods. If 
T < 2(y / 2 — 1) ps 0.83, one observes a shot noise above the Poisson noise. For r < 1 one 
has 

Pns = 4e|J| = 2P Poisson , (7.5) 

which is a doubling of the shot-noise power divided by the current with respect to the 
normal-state result. This can be interpreted as uncorrelated current pulses of 2e-charged 
particles. 

Consider next an NS junction with a disordered normal region, but with an ideal interface 
(T = 1). We may then apply the formula (3.10) for the average of a linear statistic on the 
transmission eigenvalues to eqs. (2.22) and (7.3). The result is 

(-Pns) 2 P 4 2 

(GW "32^^ (Pns) " 3 e|/| " 3 P p— ( 7 - 6 ) 

Equation (7.6) is twice the result in the normal state, but still smaller than the Poisson 
noise. Corrections to (7.6) are of lower order in N and due to quantum-interference effects 
[88]. 

Finally, consider an NS junction which contains a disordered normal region (length L, 
mean free path /) as well as a non-ideal interface. The scaling theory of section 5.2 has been 
applied to this problem in ref. [81]. Results are shown in fig. 17, where (i-Ns)/-Ppoisson is 
plotted against YL/l for various T. Note the crossover from the ballistic result (7.4) to the 
diffusive result (7.6). For a high barrier (r <C 1), the shot noise decreases from twice the 
Poisson noise to two-thirds the Poisson noise as the amount of disorder increases. 



VIII. CONCLUSION 

We have reviewed a scattering approach to phase- coherent transport accross the interface 
between a normal metal and a superconductor. For the reflectionless-tunneling phenomenon, 
the complete equivalence has been demonstrated to the non-equilibrium Green's function 
approach. (The other effects we discussed have so far mainly been treated in the scattering 
approach.) Although mathematically equivalent, the physical picture offered by the two 
approaches is quite different. We chose to focus on the scattering approach because it makes 
direct contact with the quantum interference effects studied extensively in the normal state. 
The same techniques used for weak localization and universal conductance fluctuations in 
normal conductors could be used to study the modifications by Andreev reflection in an NS 
junction. 
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FIG. 17. The shot-noise power of an NS junction (in units of -Pp i S son = 2e|/|) as a function of 
the length L (in units of l/T), for barrier transparencies T = 1,0.9,0.8,0.6,0.4,0.2 from bottom to 
top. The dashed curve gives the limiting result for F <C 1. For L = the noise power varies as a 
function of T according to eq. (7.4), between doubled shot noise ((-Pns) = 4e|/|) for a high barrier 
(r <C 1) and zero in the absence of a barrier (T = 1). For L — ► oo the noise power approaches the 
limiting value (Pns) = |e|^| for each T. (From ref. [81 J.) 

In the limit of zero voltage, zero temperature, and zero magnetic field, the transport 
properties of the NS junction are determined entirely by the transmission eigenvalues T n of 
the normal region. A scaling theory for the distribution of the T n 's then allows one to obtain 
analytical results for the mean and variance of any observable of the form A = I] n a(T n ). 
The conductance is of this form, as well as the shot-noise power. The only difference with 
the normal state is the functional form of a(T) (polynomial in the normal state, rational 
function for an NS junction), so that the general results of the scaling theory [valid for any 
function a(T)] can be applied at once. At finite V, T, or B, one needs the entire scattering 
matrix of the normal region, not just the transmission eigenvalues. This poses no difficulty 
for a numerical calculation, as we have shown in several examples. However, analytical 
progress using the scattering approach becomes cumbersome, and a diagrammatic Green's 
function calculation is more efficient. 

Note added February 1995: The theory of section 4 has been extended to non-zero 
voltage and magnetic field by P. W. Brouwer and the author (submitted to Phys. Rev. B). 



The results are 5G m (V = 0, B ^ 0) = \e 2 /h, 5G NS (V ^ 0, B = 0) = |e 2 //i, 5G NS (V ^ 



0, B ^ 0) = 0. The disagreement with the numerical simulations discussed in section 4 is 
due to an insufficiently large system size. 
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